Darkon fluid - a model for the dark sector of 

the Universe?* 

P.C. Stichel 

An der Krebskuhle 21, D-33619 Bielefeld, Germany ' 

Wojtek J. Zakrzewski 
Department of Mathematical Sciences, 
University of Durham, Durham DH1 3LE, U.K.* 

February 23, 2012 



Abstract 

We introduce darkons as fluid particles of a Galilean massless self- 
gravitating fluid. This fluid exhibits anisotropic scaling with z = |. 
The minimal gravitational coupling dynamically generates a gravita- 
tional mass density of either sign. Hence such fluid may serve as a 
model for the dark sector of the Universe. Its cosmological solutions 
give a deceleration phase for the early Universe and an acceleration 
phase for the late Universe. Will the steady flow solutions lead to a 
confining potential and so a possible model for halos? 
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1 Introduction 



Understanding the dark sector of the Universe is one of the greatest challenges 
of today's theoretical physics. Is it possible to describe the corresponding 
astrophysical observations by means of the known concepts or do we need 
some kind of 'new' physics? 

Dark matter has been introduced already in the thirties of the last century 
as the missing mass within galaxies in order to understand the galaxies' 
rotation curves and the motion of individual stars (see [1] for the history of 
dark matter). But the present model of cold dark matter (CDM) (see the 
very recent review [2]), although very successful in explaining the large scale 
experimental data, suffers from two insufficiencies: 

1) None of the proposed constituents (particles) has been observed [2], [3]. 

2) Rotation curves of dark matter dominated galaxies behave in their 
inner part in sharp contrast to the CDM-based simulations (known as the 
core-cusp problem [4]). 

Astrophysical observations of the type la supernovae and other data (see 
[5], [6] for recent reviews) suggest that the Universe is undergoing an accel- 
erated expansion. This conclusion was drawn by interpreting the data in 
the framework of the cosmological Friedmann-Lemaitre-Robertson- Walker 
(FLRW) model, which has been derived from General Relativity (GR) by as- 
suming that the Universe is homogeneous and isotropic on the largest scales 
(cf. [7]). Within this framework the cosmic acceleration is attributed to 
the so-called dark energy (DE) which is the source of repulsive gravitation 
due to its negative pressure. The simplest DE-model uses a positive cos- 
mological constant A. This concept, unified with the CDM-model, leads to 
the standard ACDM-cosmological model which is 'favoured by a number of 
observations' [5]. However, when interpreted as the energy density of the 
vacuum, the experimental value of A turns out to be a factor of about 10 121 
too small ([8], [5]). On the other hand, if we assume that DE has a dynamic 
origin (corresponding to a time dependent A) one has to change the matter 
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part of Einstein's equations by introducing e.g. a scalar field which produces 
the required negative pressure (cf. [5], [6]). Such models, however, have less 
predictive power as one can always construct a potential that gives rise to 
a given cosmological evolution. Modifying the geometric part of Einstein's 
equation is another possibility (cf [5], [6]). Such models, however, suffer from 
having to rely on an arbitrary function which is not derivable from any known 
fundamental assumption. 

In summary, so far we have no model which has been derived from funda- 
mental (conventional) physics which can explain the late-time cosmic accel- 
eration. Of course, by making this statement we assume that the accelerating 
expansion of the Universe is a real effect and not an apparent or transient 
effect due to the averaging over large scale inhomogeneities in the Universe 
(cf [9], [10]). 

Very recently, the present authors, presented a new model of the dark 
sector of the Universe, which contains no new parameters at the microscopic 
level and which is based only on some well known physical principles (Galilei 
symmetry and Einstein's equivalence principle) [11], [12]. Our model is non- 
relativistic but this in accordance with the well-known fact that the FLRW- 
model can be derived from Newtonian gravity (by heuristic arguments one 
can even include pressure [13]). 

We note that if we want to consider some new nonrelativistic particles as 
the source of the accelerating expansion of the Universe these particles must 
necessarily be massless (massive particles always lead to attractive gravita- 
tion). Nonrelativistic massless particles can exist as a dynamical realization 
of the unextended Galilei algebra. Furthermore, they live necessarily in an 
enlarged phase space [14] and possess a modified relation between energy 
and momentum (or velocity). The gravitational coupling of these 'exotic 
particles' [11] was introduced in the simplest way that satisfies Einstein's 
equivalence relation. This coupling, for a self-gravitating fluid, leads to a 
dynamically generated active gravitational mass density of either sign which 
then becomes a source of the corresponding gravitational field. This fact 
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leads to the possibility of using our model as a basis of further models which 
can then constructed to explain the observed accelerating expansion of the 
Universe and also of the CDM-attributed halos which confine the galaxies. 

The equivalence principle, i.e. the local equivalence between gravita- 
tion and acceleration is normally formulated mathematically as an invari- 
ance w.r.t. arbitrary time-dependent translations (gauge transformations). 
By fixing this gauge we can reduce the phase space and remove the unphysi- 
cal variables. The Lagrangian of the resultant model does not possess a free 
particle limit and hence the particles it describes, darkons, exist only as fluid 
particles of a self-gravitating fluid. The darkon fluid presents a dynamical 
realisation of the zero-mass Galilean algebra extended by anisotropic dilation 
symmetry with the dynamical exponent z — |. 

The paper is organised as follows. In Section 2 we introduce our Galilean 
massless particles and discuss their coupling to gravity. In Section 3 we 
describe the generalization of this particle picture to a self-gravitating fluid. 
In Section 4 we consider the cosmological solutions and in Section 5 we state 
the main open problems. We conclude in Section 6 with some final remarks. 

2 Galilean massless particles and their cou- 
pling to gravity. 

Our main task involves finding the minimal dynamical realization of the 
unextended (m = 0) Galilean algebra in d = 3. To do this we consider the 
first order Lagrangian in a 12-dim phase space 

L = Piii + qiVi -ptVi (1) 

leading to the equations of motion (EOM) 

Xi = y%, Pi = 0, qi = -pi, and y 4 = 0. (2) 

Here x iy resp. yi are the components of the position, resp. velocity, 
of a particle and p iy resp. q iy are the corresponding canonical momenta. 
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They satisfy, in the framework of the equivalent Hamiltonian dynamics, the 
Poisson-bracket algebra (PBs): 



[xi, Pj] = Sij and [y h qj] = 



(3) 



with the Hamiltonian H given by 



H 



PiVi- 



(4) 



We note that 

• The relation between energy H and the momenta (velocities) is non- 
standard, 

• Momenta Pi and velocities yi are dynamically independent, 

• The qi have no counterpart in standard particle mechanics. 
By Noether's theorem we obtain from (1) the boost generator 



which, according to (3) has vanishing PBs with the generator of the spatial 
translations Pj — pf 



So (6) clearly shows that the Lagrangian (1) describes a dynamical re- 
alization of the unextended (m = 0) Galilei algebra Go- To prove that this 
realization of Go is the minimal one we refer to [14]. We note that Go may 
be enlarged by dilations D with arbitrary dynamical exponent z. 

D is given by [14]: 



To couple our massless particles to the gravitational field strength g^x, t) 
we have to be consistent with Einstein's equivalence principle: Locally the 
gravitational field strength is equivalent to an accelerating frame. 



Ki = Pit + qi 



(5) 



[K it Pj] = 0. 



(6) 
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So, the EOM 

Xi = yt = 

for a free particle have to be replaced, in the presence of gravity, by 



(8) 



Xi(t) = gi(x{t),t). (9) 

This EOM can be realised if we add to L an interaction part (minimal 
coupling) 

Lm = -qi9i- (10) 

Then the EOM 

ft = -Pi = (11) 

gets replaced by 

ft = -Pi = qkdig k . (12) 

Note: Einstein's equivalence principle still holds if we add in (10) further 
terms containing spatial derivatives of gi ie by replacing 

<7« 9i — Ki Agi + higher order terms (13) 

where A is the Laplacian. 



3 A self-gravitating (darkon) fluid 

In the Lagrangian formulation we generalise the one-particle phase space 
coordinates Ai £ (xi,pi, yi, qi) to the continuum labelled by £ G R 3 (comoving 
coordinates) 

The Lagrangian for the self gravitating fluid then becomes (after elimination 
of pi and yi) 

L = -n J d 3 £ (qi&^Xi&t) + qi{lt) gi {x(£,t),tfj + L+, (15) 
where, as usual, 
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^ J d 3 x(g t (x,t)) 2 . (16) 



87TG 

G is Newton's gravitational constant and uq is the (constant) particle 
density in £-space. Euler-Lagrange EOM: 

Varying g« and Xi leads to the continuum generalisation of the previous 
equations (9) and (12) for Xj and (ft. 

Now we have two possibilities to go further: 

In the first case we take <ft = —d^ and vary <fi getting 

A0 = 4nGd k (nq k ) (17) 



with 



and 



n(x,t) = n J 1 (particle density) (18) 



J = det y X Q^ J (Jacobian), (19) 

where £(x,t) denotes the function inverse to x(£,t). 

The right hand side of the Poisson equation (17) describes a dynamically 
generated active gravitational mass density which may be of either sign with 

+ sign leading to an attractive gravitation 

- sign leading to a repulsive gravitation. 

This promotes the self-gravitating fluid to a possible candidate for the 
dark sector of the Universe. 

In this case our EOM are invariant w.r.t. arbitrary time-dependent trans- 
lations ( gauge transformations) 

Xi -> Xi + Oi(t), (20) 

where transforms as 

<j>{x,t) -> <p\x,t) = <j)(x,t) - a^Xi (21) 
and q,i and n remain invariant. 
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In the second case we vary g { instead of and this leads to a linear relation 
between qi and ^ 

9i (x,t) = -4irGn(x,t) qi (£(x,t)t). (22) 
Note that the expression for ^ can be rewritten as 

Qidt) = (^)m,t),t) (23) 

The expression (23) allows to eliminate from the EOM the unphysical 
phase space variables qi in favour of the gravitational field strength gi. 

This removes the unphysical degrees of freedom. However, this expres- 
sion for qi (a particular solution of the Poisson equation) destroys the gauge 
symmetry ie it fixes the gauge. 

Remark: The gauge symmetry can be restored by considering a curl-free 
velocity field 

u = Vw 

and treating u as the variable to be varied [12]. 

Inserting q^ from (23) into the previous EOM we get two coupled EOM 
for x(£,t) and g(x(£,t),t) 

d 2 

^2 = 9M£,t)t) (24) 



(initial condition x(£ 0) = |) and 

^ (f ) (x(it),t) = (^f) m, t ),t). (25) 



The last equation can be integrated once. To do this we note that 

a ^^_dxkd_[9k\ , 9fi x 

1 ~ 3£ n d&dtKn) [ ] 

is conserved (ie j-fli = 0). 

This conservation law gives the required once-integrated form of (24) 

d (9i\ g k ~ . a d£ k 

Jt U = n 9lXk ~ 9k dx~' (27) 
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which is the transport equation for the field strength g (cf. [15] for the 
standard Newtonian case). 

Note that 

Q = J d 3 t6 k (£)a k (£) (28) 
generates infinitesimal relabeling transformations 

£ ->• £ + a(£) with V f a = 0. (29) 

Remark: 

The EOM (24), (27) do not contain the gravitational constant G. By 
inserting (23) into (15) we see that the Lagrangian contains G only as a 
common factor. The Lagrangian does not split into free and interacting parts. 
Thus the fluid particle (called darkons) do not exist as free particles but are 
only elements of a self-gravitating fluid. This comes from the reduction of 
the phase space (expression of in terms of g t by (23)). 

For brevity we dispense with the Eulerian formulation of the darkon fluid 
dynamics. Therefore we will only summarize the resulting space-time sym- 
metries [12], [14]: 

The symmetry algebra becomes the expansion-less conformal Galilei al- 
gebra with z — |. For free particles dilation symmetry can have any z, the 
minimal gravitational coupling fixes z to z — | (higher-order derivative terms 
(13) change this). 

4 Cosmological solutions 

We consider the darkon fluid as a model for the dark sector of the Universe. 
Hence we look for solutions of the fluid EOM satisfying the cosmological 
principles (the Universe being isotropic and homogeneous on large scales). 
This gives in the Lagrangian formulation 

n(x,t) = n{t), Xi&t) = a(t)&, (30) 
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where a(t) is the cosmic scale factor, 
and 

9i(x,t) = -An n(t) g(a(t)) Xi (31) 

Note that the derivative terms in the gravitational coupling would not con- 
tribute!. 

If we next put all this into the EOM we obtain: 

1. From the definition of n(x,t) 

»(*) - ^ P2) 

2. From the EOM (24) and using (32) 

a(t) = -^g(a), (33) 

Qj 

i.e. a Friedmann-like equation 

3. From the EOM for gi and the definition of 

9«t)) = -J^, P =const, (34) 
where (3 is related to by 

^(1) = ATT^i- 

The last two equations, after integration, give us a cubic equation for g{a) 

[11] 

g(g 2 + d) + Co (l - -) = 0, (35) 
with Cj integration constants and 

at = ^Co- (36) 
Next we take Q > and use a(t) as a measure of time (given that now 

we have a(t) > (expanding Universe)) 

Then we conclude from (35) and (33) that 

1. For a < a t we have g(a) > and so a < - ie the deceleration phase 
of the early Universe. 

2. For a > a t we have g(a) < and so a > and we are in the accelaration 
phase of the late Universe. 
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Notes 



• a t defines the point of transition between the phases. All this is con- 
sistent with cosmo logical observations, (cf. [11]) 

• We need C{ > for these results (this has to be explained by some 
physical arguments) 

• If Cj = we reproduce the scale invariant solution valid at small t 
(then a(t) ~ £§). 

5 Open problems 

The main open problems are the following: 

• Steady flow solutions. Question: Can a steady darkon fluid flow be a 
model for halos? ie provide a gravitational potential for the formation 
of stable structures (galaxies or clusters of galaxies)? 

For that we need to find steady solutions of our EOM leading to a 
potential <fi which confines massive objects. Our preliminary results 
for the spherically symmetric case demonstrate the existence of such a 
confining potential at least for small distances from the centre. But to 
get a quantitative agreement with observations we, probably, have to 
include spatial derivative terms in the gravitational coupling (cf. eqs 
(10) and (13)). 

• Embedding of the model into a larger one which possesses a general 
relativisitic limit (as z — | our model cannot be a non-relativistic limit 
of a realtivistic one (z = 1)). 

6 Final remarks 

Our main results are: 
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• Parameter free Lagrangian for a self-gravitating system of Galilean 
massless particles (darkon fluid) 

• Dynamical generation of an active gravitational mass density of either 
sign. 

• Explanation of the deceleration phase of the early Universe and of the 
acceleration phase of the late one. 
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